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Abstract: 

The periodic orbits of the strongly chaotic cardioid billiard are studied by introducing a binary 
symbolic dynamics. The corresponding partition is mapped to a topological well-ordered symbol 
plane. In the symbol plane the pruning front is obtained from orbits running either into or 
through the cusp. We show that all periodic orbits correspond to maxima of the Lagrangian 
and give a complete list up to code length 15. The symmetry reduction is done on the level of 
the symbol sequences and the periodic orbits are classified using symmetry lines. We show that 
there exists an infinite number of families of periodic orbits accumulating in length and that 
all other families of geometrically short periodic orbits eventually get pruned. All these orbits 
are related to finite orbits starting and ending in the cusp. We obtain an analytical estimate of 
the Kolmogorov- Sinai entropy and find good agreement with the numerically calculated value 
and the one obtained by averaging periodic orbits. Furthermore the statistical properties of 
periodic orbits are investigated. 
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1 Introduction 



A key step towards an understanding of the behaviour of a dynamical system is achieved by 
finding a symbohc dynamics. By means of the symbohc dynamics trajectories can be labeled 
by doubly infinite sequences built from a (not necessarily) finite alphabet. Periodic orbits are 
represented by periodic sequences and can be systematically searched for once the coding is 
known. The knowledge of a complete set of a large number of periodic orbits up to a given 
geometric length is necessary for the application of Gutzwillers periodic orbit theory [1], which 
relates the quantum mechanical density of states of the quantized billiard system to a sum over 
classical periodic orbits. If the symbolic description is not complete there are sequences which 
do not correspond to physical orbits. To deal with this more complicated case the idea of a 
pruning front in the symbol plane was introduced in [2]. These methods have been successfully 
applied to a number of systems. There are, however, no systematic methods to find a symbolic 
description and the corresponding pruning mechanisms. 

In the class of billiards inside simply connected domains of the Euclidean plane ergodic examples 
typically have either families of orbits accumulating in length, singularities in the boundary or 
non-isolated parabolic families. In the case of the cardioid billiard, which has been rigorously 
proven to be strongly chaotic, i.e. it is ergodic, mixing, a A^-system and even a Bernoulli 
system [3, 4, 5, 6], we have accumulating families and one singularity. The relation of the two 
in the cardioid billiard is quite interesting and a thorough understanding of their effects is a 
prerequisite for the semiclassical quantization of this system. 

There already are several results for both, the classical and quantum mechanical billiard. The 
cardioid is the limiting case of a family of billiards introduced by Robnik [7], see also [8, 9] 
and references therein. The statistical properties of the eigenvalues of the quantized cardioid 
billiard were studied in detail in [10, 11]; see also [12], where the focus is on diffraction effects. 
A lot of work has been done on Robnik's family but the classical mechanics of the cardioid has 
not been analyzed in depth. It is this gap we want to fill in with this work, keeping in mind, 
however, the application to quantum mechanics. 

The paper is organized as follows. In sec. 2 the cardioid billiard and the billiard map are defined. 
We show that products of linearized maps always have a positive trace. Therefore all orbits are 
maxima of the Lagrangian. Subsequently a discussion of the symmetries of the billiard map 
is given. In sec. 3 the symbolic dynamics is defined and the corresponding partition of the 
Poincare section is illustrated. The initial partition is given by the discontinuity of the map. 
The pruning of code words is discussed in the symbol plane using the pruning front, which 
is also obtained from the discontinuity of the map. In sec. 4 the periodic orbits are classified 
according to their symmetry by using the symmetry lines of the desymmetrized billiard. We 
give a list of the number of periodic orbits in each symmetry class up to code length 15. Families 
of periodic orbits with short geometric length and their relation to cusp orbits are investigated 
next. It is shown that most of them eventually get pruned. However, there remains an infinite 
number of families accumulating in length. Concerning the number of cusp orbits we show that 
they are more frequent than periodic orbits. In sec. 5 we obtain an analytical estimate for the 
Kolmogorov- Sinai entropy and find good agreement with the value obtained from numerical 
methods. The average length and the KS entropy are calculated using the periodic orbits. 
Finally, we focus on the statistics of periodic orbits. 
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2 The cardioid billiard 



A billiard inside a two dimensional Euclidean domain is given by the free motion of a point 
particle inside with elastic reflections at the boundary 90, i.e. the angle of incidence equals 
the angle of reflection. The cardioid billiard is the limiting case of a family of billiards first 
studied by Robnik [7]. Their boundary in polar coordinates (/?, (f>) is given by 



p{(f>) = f + e cos ^, 



G [-7r,7r]. 



We restrict our attention to the cardioid (see fig. f) which is obtained for e = f , or implicitly 
by 



0, 



F{x,y) 

where (x, y) = r[(f>) = {p{<f>) cos ^, p[(f>) sin (f>). At (f> = zLtt the cardioid has a singularity located 
at the origin r(ib7r) = (0,0). 
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Figure f : The full and the desymmetrized cardioid billiard. 
From the above definition (f ) one easily derives the curvature 



the unit tangent vector 



r = (r.,r,) 

ds 
d<t> 



4cos((^/2)' 
sin(3(^/2),cos(3(^/2)). 



(3) 



(4) 



and the differential of the arc length = 2cos(^/2). Thus the arc length s is related to (f) by 



s = 4sin((^/2) . 

The area of the cardioid is |0| = 37r/2 and its circumference is \dft\ = 8. 



(5) 



2.1 Poincare map 

We now derive the Poincare map from bounce to bounce in Birkhoff coordinates ^ = (s,p) where 
p is the component of the velocity in the direction of T right after the refiection. The Cartesian 



3 



components of the unit velocity u of a billiard ball starting on dft at r[(f>) is determined by 
the angle /3 G [pz/2,7r/2] measured with respect to the normal N = (—Ty^Tx) pointing inward; 
/3 = — 7r/2,0,7r/2 specifies the direction — T, T respectively. The velocity in the T^N 
coordinate system is denoted by (p, n) = (sin /3, cos so that we obtain v = (— cos(/3 — 
3^/2),sin(/3 — 3^/2)). The rightmost point of dft corresponds to arc length s = 0, so that 
s extends from —4 to 4 while p G [—1,1]. We now construct the Poincare map P from the 
rectangle V = [—4,4] X [—1, 1] onto itself by calculating the image of a point ^ = G V. 

We seek the first intersection of the ray starting at r = (x, y) in direction v with the boundary 
curve dft 

F{x + tv,,y + tvy) = 0, F{x,y) = 0. (6) 

One root of this fourth order polynomial in t is always zero, because of F[x^y) = 0. Thus we 
obtain a polynomial Q{t) of degree 3 for t 

Q{t) = t^- (^4cos(-/3 + ^) + 2cos(/3 + ^)^t2 (7) 
+ (7 + 8 cos{<p) + 4 cos(-2/3 + <p) + cos(2/3 + <p) + i cos(2/3)) t/2 
-8cos(/3)cos3(-^) . 

The smallest zero t* > is the physically meaningful solution (only for the sliding motion we 
have t* = 0). To complete the description of the Poincare map we first calculate the new point 
r' = (x', y') = (x + t*Vx^ y + t*Vy) on dft and determine 

(f>' = arctan(j/'/a;') 

p' = sin{^') = {T',v') = {T',v) = sin{3{<P'-<P)/2 + ^) . (8) 

The complete map P : ^ = (s^p) ^' = {s\p') is now given by the composition of the above 
steps. P is area preserving because s and p are canonically conjugate. For convenience we 
sometimes use (f> instead of s, but without mentioning we assume it to be expressed in terms 
of s. 

The two curves 

S± = {CeV\s = ±i}, (9) 

corresponds to orbits which start in the singularity (f> = zLtt, i.e. on the right or left boundary 
of V. For the image of points ^ = (s,p) G S± under P one has (f>' = aTctan[vy/vx) because of 
r = (0, 0). Thus one has (f>' = ib7r/2 — /3 and using equation (8) p' = — sin(ib7r/4 — /3/2). The 
two image curves join at the origin of V. For reasons that become clear later on, we denote 
them by 

Tg' = {UV\p=-s/i, ±s>0} , (10) 

such that 

P{S±) = Tl\ (11) 

Fach of the curves iS+ and S- has a fixed point (4,-1) respectively (—4, 1), which corresponds 
to Vx = —1. The physical motion starting at these fixed points is a sliding motion along the 
boundary, either counterclockwise or clockwise. 

Note that although r(— tt) = r(7r) = (0,0) we take them as different points equipped with 
their (different) tangent vectors. In a differentiable point of dft one can start in directions 
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/3 G [— 7r/2, 7r/2]. In singular points this interval can be different; in our case at r = (0,0) we 
can go in any direction. Points starting with Vy > are attached to ^ = tt, the ones with Vy < 
belong to ^ = — tt. The unique point with r = (0,0), v = (1,0) in phase space can be assigned 
either of the two coordinates (4, 1) or (—4, —I) in V. This is only a coordinate singularity and 
P correctly maps both points onto the same point ^ = (0,0), where the images of S± meet. 
Strictly speaking the map P is defined on a rectangle with two opposite corners identified. 

Two other special lines in V are its upper and lower boundary 

J'± = {UV\p = ±l,s^±i}. (12) 

Note that J-'± are half open intervals. J-'± defines starting points outside of the singularity with 
a velocity parallel to T, i.e. /3 = ib7r/2, p = ±1. All the points of J-'± are fixed points of the 
map, although physically they correspond to the above mentioned sliding motion. In a physical 
sense all of them are images of the two points (±4, =pl) where JF, S and intersect (omitting 
the index it refers to both lines). 

Reversing the velocity on we define the line 

T = {UV\p = s/i}, (13) 

which is the set of initial conditions that will right away hit the singularity. Therefore is 
the mirror image of F with respect to the s-axis. F is of utmost importance, because this line 
turns out to be 1) the discontinuity of the map, 2) the boundary between our primary symbol 
regions and 3) the origin of the pruning front. F separates two regions A and B in V (see 
fig. 2), where A is below F and B above. We consider A and B as open sets, i.e. without the 
lines S± and F. 

1 

P 


s_ 
-1 

-4 -2 F_ 2 s 4 

Figure 2: The regions A and B separated by F in the Poincare surface of section togeth- 
er with the fix lines J-'± and the singularity lines S±. The shaded regions correspond 
to four intersections of the line r -\- tv with the boundary 90, i.e. three real solutions 
of the polynomial Q{t). In the white regions there always is only one real solution. 
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Figure 3: Demonstration of the kneading property of the biUiard map. The left picture 
shows the regions A and B which are mapped by P to give the picture on the right 
hand side. P "kneads" the triangle B by sliding the point 4 to the lower right corner, 
while in turn the point 1 is taken to the middle. The remaining points 2 and 3 stay 
fixed. Note that the former center 4 becomes a corner and vice versa for 1. In a similar 
way A is deformed under P, by taking 5 to the upper left corner, and moving 6 to the 
middle. The mapping P thus contracts in the direction of F, and expands along F~^. 

In order to understand the "kneading properties" of P we need to know the behaviour of P near 
F. Since P is discontinuous on F there are two different continuations of P on F. Approaching 
F from a region with negative discriminant of Q{t) (shaded regions in fig. 2) the limit orbits 
hit the singularity and refiect at the horizontal limiting tangent. We denote the corresponding 
limit map by P^ and find (it denotes the sign of s) 

P,(5,5/4) = (±4,±I-5/2) . (14) 

In the other case we ignore the encounter with the cusp and find 

Pj(3,3/4) = (-(±4-3),(±4-3)/4). (15) 

In fig. 3 the kneading property is illustrated by mapping the regions A and B. Here the lines 14 
and 56 map according to Pj, while 85 and 43 map by P^. The main tendency is an expansion 
along the direction of F~^ and a contraction along the direction of F. 

Since the map is discontinuous at F, the image of F under the Poincare map is not well 
defined. Moreover, there is no unique tangent vector in iS, because the boundary curve is not 
differentiable in this point. Note that we can assign the coordinate s' = ±4, but we cannot 
specify p' for G F. One might be tempted to define the image of F according to its limit under 
Pj, because also the corresponding tangent vector in S± can be defined by an appropriate one 
sided differentiation. This is, however, misleading as we will see in the discussion of finite orbits 
that start and end in S. 

If we consider the desymmetrized billiard, see fig. I, the corresponding billiard map P defined 
on V = = (s,p) I s G [0,4], p G [—1, 1]}, is obtained by first using P for a given ^ ^ V; if 
^' ^ V we have ^' = i^', otherwise ^' = (— s', —p'). 

In the numerical calculation of P there is one difficulty (besides hitting S exactly) if we are 
close to F, especially around ^ = (0,0). Firstly the discriminant of Q{t) changes quite rapidly 
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and is very close to because we are close to a double root ol Q{t). Secondly at the new point 
of dft we have very small r and the calculation of s' is susceptible to round-off errors. The 
reason is the unbounded curvature close to the singularity of dft. 



2.2 Linearized map 

For an arbitrary billiard the linearized Poincare map from i^i to ^2 can be expressed as (see, 
e.g. [13]) 

= ( ^^0' I ) ( ? ) ( 0^ -lln, ? ) ( ''o IM ) ■ ^^^^ 

Here Hi and rii denote the curvature and normal component, respectively, and / is the geometric 
length between the two reflections. Lower indices indicate that the corresponding quantity is to 
be evaluated at Si. Double indices for quantities depending on both points are usually omitted. 
For a periodic orbit this reduces to the more familiar form for the monodromy matrix M (see, 
e.g. [14]) 

^^^n(; 'r )(2;;„. -0- <"> 

We prefer to use (16) because it has a special form which allows us to show that every peri- 
odic orbit (except on JF) has positive trace. The argument is along the lines of Wojtkowski's 
pioneering work [3], but we will use a slightly different form due to Wittek [15], who applied 
it to the wedge billiard. The key point is that -DP21 always has the following checker board 
structure (excluding JF) 

<-) 

Note that the product of two checker board matrices is again a checker board matrix. Let us 
rewrite (16) 

DP -i^l''' ° U ^'-^ -I \(n, \ 

Since rii > and / > (excluding J-'±) we need to show ki > 1 and 

Because of Hi > we have ki > and therefore the latter inequality implies ki > 1. Fvaluating 
the trace with these inequalities we find that Tr _DP2i ^ 2, where equality holds for ki = k2 = 1. 
Starting from the explicit description of the Poincare map in the previous section it seems quite 
hard to obtain this statement because the solution of a cubic equation is involved. Therefore 
we now look at the generating function of our map (see [16] for a review), which is just the 
length between successive refiections at r(^i) and r(^2)- Note that in order to obtain the area 
preserving map P in coordinates = (s,p) we should parameterize by Si and S2- However, the 
calculations are more conveniently done with /(^i, ^2)- Denote by L the vector joining the two 
refiection points, 

L{<f>r,<f>2) = r{<f>2)-r{<f>r) (21) 
^1,(^2) = (22) 



2\lsin j/(cos2 J/ -|- 2 cos x cos J/ -|- 1) (23) 
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where x = [(pi + ^2)/2 and y = [(pi — ^2)/2. The unit velocity is given by v 
obtain 

ni = {Ni,v) n2 = -{N2,v), 



= L/l and we 
(24) 



cos^ y + 2 cos X cos y + I 



(25) 



(26) 



2 cos^ J/ + 3 cos x cos J/ + 1 lb sin x sin y ' 



where + and — corresponds to ki and A;2, respectively. Finally 



1 1 

ki k2 



2 2 cos^ J/ + 3 cos x cos j/ + 1 



(27) 



3 cos^ J/ + 2 cos X cos j/ + 1 



1 sin^ 
1 ^ 



3 cos^ J/ + 2 cos x cos y + I 



(28) 




(29) 



proves the inequality (20). The geometric origin of this relation is the convex scattering property 
[3] of the cardioid, i.e. d^n/ds^ > 0. There are a some very important consequences: 

• The eigenvalues of the monodromy matrix of periodic orbits are always positive, i.e. all 
periodic orbits are direct hyperbolic. We conclude that periodic orbits are maxima of the 
Lagrangian C = /(^i-i, ^i) [17, 16]. Thus the numerical search for periodic orbits is 
very much simplified, because we do not have to find saddle points of £, which would 
correspond to inverse hyperbolic orbits (see [18] for related results). 

• The maximum number of conjugate points along a periodic orbit is given by the number 
of refiections. The reason for this is contained in the optical interpretation of (20) already 
described by Wojtkowski [3]: there is enough time between successive refiections in order 
for a conjugate point to occur. Since for a free motion there cannot be more than one 
conjugate point, the above statement follows. 

• We can obtain an analytical estimate (from below) of the maximum Lyapunov exponent 
and therefore also for the Kolmogorov- Sinai entropy by a theorem from [3], see sec. 5. 

Note that the first statement only holds in the non symmetry reduced system. This is one 
of the reasons why we think that for the cardioid it is worthwhile to study the non-reduced 
system. Furthermore notice that the convex scattering property, which is the basis for all of the 
above, does not hold for any other member of the family of billiards (1), which are therefore 
much more difficult. 

2.3 Symmetries 

The time reversal symmetry of a billiard combined with the spatial symmetry of the cardioid 
gives us a number of symmetry classes of orbits. We will not pass to a desymmetrized billiard 
map but instead do the symmetry reduction on the level of the symbolic dynamics. Here we 
discuss the manifestation of the symmetries in the map. In the following sections this will be 
translated into symbol plane. 
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Figure 4: The symmetry lines in the Poincare section V . Tq corresponds to orbits 
starting at right angles on the boundary, Xq contains all orbits starting in s = 0, and 
corresponds to orbits intersecting the at right angles. 

The standard time reversal symmetry is given by 

T ■.{s,p)^{s,-p) . (30) 
As usual the T-symmetry allows for a simple expression of the inverse map as = TPT . 



The cardioid billiard also has reflection symmetry with respect to the x-axis 

X : {-s,p). 
X and T are involutions, i.e. = id, X"^ = id and det T = — f , det X = —1. 



(31) 



Furthermore we have P 
involutions 



XP^'X and P- 



TP^T^ and thus we define two families of 



T 

X,n 



rp2 

n 

xl 



id 
id 



det Tn 
det X,n 



n = ±f,±2... (32) 
n = ±f,±2... (33) 



The fixed point sets of involutive symmetries, the so-called symmetry lines, are useful in finding 
symmetric periodic orbits [f9] and in their classification [20]. We define 



= {C\TnC = C} 

The set Tq contains all the orbits starting at right angles on the boundary, i.e. 

To = {UV\p = 0}. 
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(34) 
(35) 



(36) 



Xq contains all orbits starting at the rightmost point of the cardioid with arbitrary angle 



Xo = {CeV\s = 0}. (37) 
All iterates of the symmetry lines can be obtained by iterating Tq^ Tl, Xq and A*!, because [19] 

^2n — -T -!0 -'2)1+1 — -T 

y pra y y nn -y v*^"/ 

A.2n — -T ^0 ^2n+l — r A-i . 

The equations PT^^ = ^ for Ti do not have a solution in billiards without potential and therefore 

also T2n+1 = 0. 

The symmetry line = {(^ | PX^ = ^ : (— = corresponds to orbits that intersect 

the at right angles. Therefore we introduce the suggestive notation 

'^2n = '^2n (39) 
= ^2n+l (40) 

related to the geometric form of the corresponding periodic orbits. is obtained by using (8) 
which yields the condition sin(3^ + /3) = sin(/3) and finally 

4+ = {e|p = -cos(3(^/2),(^G [0,7r/2]} (41) 
4~ = {e|p = cos(3(^/2),(^G [-7r/2,0]} , (42) 

because for \(f>\ > 7r/2 we cannot cross the x-axis. 

If we consider the family of billiards (1) the set S does become a symmetry line. In the spirit 
of this notation we should call it X^. Moreover, the image respectively preimage of F = A!^2 
respectively = is well defined in this case. The interpretation of the line of discontinuity 
r as a symmetry line leads to a nice interpretation of finite orbits starting and ending in S . 

At intersections X^ H X^^ X^ H 7^ and 7^ fl 7^ of symmetry lines there are periodic orbits. 
Intersections of symmetry lines of the same type X oxT are periodic orbits with (not necessarily 
primitive) period I m — n I [19]. If A'^ intersects 7^ we instead have period |2n — 2m |. This is easily 
seen in the following way: We have P'^T^ = and P'^X( = thus one has p-^-^TX( = 
which in turn implies p2(ra-m)^ _ After symmetry reduction we always find period |m — n\. 

Let us now study the the intersections of symmetry lines visible in fig. 4. Already the basic 
lines Tq and Xq intersect at ^ = (0,0), but the corresponding orbit is not periodic but instead 
a finite orbit running along the x-axis. The point ^ = (0, 0) is also an intersection of F and 
F~^, i.e. the corresponding orbit starts on S, is mapped to F~^ and F and back into S (with 
undefined p'). 

The periodic orbit with the shortest period (besides all the parabolic fixed points on J-'±) is 
given by the intersection of Tq and X^^ at (±2,0) respectively <f> = ib7r/3. There is one more 
intersection in fig. 4: X^ has endpoints on F. The (well defined) preimage of this point on F i 



IS 



on F~^, so that this orbit has two refiections besides the point in the singularity. A discussion 
of symmetric orbits with higher period will be postponed until we have the symbolic dynamics 
at hand. 
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Figure 5: Examples for the determination of the symbols A and B in position space by 
the following rule: Connect the singularity S with the current point P. Now determine 
if the velocity vector is inside the sector formed by PS and the oriented tangent vector. 
If the velocity vector is inside the sector then B is assigned, otherwise A. In the left 
example, the symbol B is associated, whereas in the right example, the corresponding 
symbol is A. 

3 Symbolic Dynamics 

The discussion of the Poincare map and its discontinuities proposes a natural choice for the 
initial partition of V: the regions A and B as separated by F. For a given velocity on dft 
one can easily read off its symbol as illustrated in fig. 5. Transforming this description in the 
Poincare section to the one in configuration space means to consider two consecutive points 
(f> and (f>' on the boundary. If ^' > ^ the letter B is assigned for ^, if ^' < ^ we obtain A. 
Furthermore we exclude the cases cf) = cf)\ cf) = zLtt and cf)' = zLtt in order to ensure that the 
regions A and B correspond to open sets in V . 

The construction of the partition is as follows: If we superimpose F and F~^ we obtain four 
cells in V labeled A. A, A.B, B.B, and B.A which are shown in fig. 6. The forward image of 
these cells generates "past" stripes AA., i?A., AB. and BB. basically along the direction of 
F~^ ordered by increasing p. The operation of P on a symbol sequence (word) just shifts the 
dot to the right. The sequence AB. for a strip means that its preimage is in .B and moreover 
in that part of .B whose preimage is in .A. 

The backward images of the initial cells A. A, A.i?, B.A, and B.B give stripes elongated in the 
direction of F labeled .AA, .AB, .BA and .BB, which are just the images of the past stripes 
under T . The new words .AA, .AB, .BA and .BB tell about the future of a given strip. The 
intersection of these two sets of stripes generates a partition of V into 2'^ cells which is also 
displayed in fig. 6; each cell is uniquely labeled by four symbols. Iterating this process gives a 
finer and finer subdivision, where the hope is that the cell size eventually goes to zero so that 
every point in V can be uniquely labeled by a doubly infinite sequence. We cannot prove that 
this process generates a partition, but we have good numerical evidence that it does. 



11 




< < <■ < CQ tti ffl ffl' 
<<fflffl<<fflffl 




Figure 6: The division of V by preimages of V and images of into 2^ cells with 2 
symbols, 2'^ cells with 4 symbols and 2^ cells with 6 symbols. On the right hand side the 
corresponding division of the symbol plane is shown. Every line has 1, 3 respectively 
7 intersections with other lines. 
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Figure 7: The division of V into 2^° — 4 cells with 10 symbols. At this iteration level 
four lines have less than 31 intersections, indicating the onset of pruning. 

3.1 Symbol plane 

For any word a.u = . . . s_2S_i .S1S2 • • • a point ((5, 7) G [0, 1] X [0, 1] in the symbol plane for the 
"future" -coordinate 7 is calculated by 

00 

7 = E^»-2"'' . (43) 

8 = 1 

and for the "past" -coordinate by 

00 

8 = Y,s_,2-^ , (44) 

8 = 1 

where in the context of numerical interpretation Si is zero (or one) for symbol A (or B). 

A quite surprising observation is that the ordering of stripes in the Poincare section corresponds 
to the ordering of words ((5, 7) in the symbol plane. Thus the symbolic dynamics is already 
well-ordered [2, 21]. The dynamics in the symbol plane is a shift on the symbol sequences, i.e. 
it maps according to the baker map (see, e.g. [22]). 

Fig. 7 shows the onset of pruning. For symbol length 8 we have 2® cells, i.e. the 16 past and 
16 future stripes intersect pairwise. At symbol length 10 this is not true any more. The future 
stripe .AB'^, which is the first one below F does not intersect the two outermost past stripes 
i?^., so that there is no cell e.g. with label B^ .AB'^ (see fig. 8 for a magnification of the relevant 
region in V). In the magnification it is also visible, that the future stripe .BA^ above F still 
intersects the past stripe i?^., such that the cell B^ .BA^ exists. 
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Figure 8: Magnification of fig. 7. Tlie upper riglit stripe is B^., tlie stripe above F is 
.BA^, and .Ai?'* is tlie one below. Since the stripes B^ . and .AB'^ do not intersect, the 
sequence B^.AB"^ is forbidden, whereas .BA"^ and B^ . still intersect, such that B^ .BA^ 
is allowed. 




0.5 5 1 



Figure 9: 10^ iterations of a single orbit shown in the symbol plane truncated to word 
length 2x20. 

To find out systematically which cells are forbidden we plot the truncated symbolic past and 
future for points from a generic trajectory in the symbol plane (see, e.g. [21]). The result is 
shown in fig. 9. One observes that the picture is symmetric under refiections with respect to 
the diagonals. This is caused by the underlying symmetries T, X and TX: For a word a.u 
the symmetry operations are realized by the operations of reading backwards to. a and taking 
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the complement a.ib. T is realized by cu.a, and X is represented by to. a alone. Thus TX 

is represented by a.ib. Using eqs. (43) and (44) one can easily calculate the action of these 

operations in the symbol plane. For a given word a.u with coordinates ((5, 7) in the symbol 
plane we have 

Ta.u = u.a T{S,-f) = (1 - 7, 1 - (5) 

Xa.uj = u.a ^('^,7) = In^S) (45) 

TXa.u = a.u TX{6rf) = (1 - (5, 1 - 7) . 

Similarly we now can define the basic symmetry lines for the symbol plane. is given by the 
diagonal 8 = ^ and Tq is the other diagonal (5 + 7 = 1. A detailed discussion of the symmetries 
in connection with periodic orbits will be done in sec. 4. 

3.2 Pruning front 

Returning to the onset of pruning at word length 10 close to V one expects from the geometry 
of the billiard that the origin of pruning is the singularity of the cardioid. Therefore, we now 
study the set of symbol sequences of orbits that almost hit the singularity. Two possibilities 
arise: on the one hand the orbit may just miss the singularity, or on the other hand, it may be 
refiected very close to the singularity. The limiting cases of these types of orbits generate the 
pruning fronts [2, 21] between allowed and forbidden sequences, see fig. 10. 

We think that this pruning front fulfills the conjecture stated in [2] that the region enclosed 
by the front and the (5-axis specifies the primary pruned region in the symbol plane. All the 
other forbidden cells visible in fig. 9 are related via the symmetry operations or they are images 
or preimages of these regions. We cannot prove that the pruning front separates allowed and 
forbidden orbits in the symbol plane, but we have found no counterexample. 

The two pruning fronts are obtained from the two possible limiting maps Pj and on V. 
The recipe is to start with G F and to use either Pj or P^ to map across the discontinuity. 
The corresponding symbols can be read off from figures 2 and 3; if we start for example close 
to the line 56 shown in fig. 3, the current symbol is A and the application of Pj carries us 
into B. In the case of i^' G <S we assign A if s' = 4 and B ii s' = —4 for the refiection in the 
singularity. After this crucial step all further forward images under P are well defined, similarly 
the preimages of G F. Therefore an infinite symbolic past and future can be assigned to 
and the corresponding point in the symbol plane is part of the /- or s-pruning front, depending 
on the initial mapping step. From the two pruning fronts in fig. 10 the left with 8 < 0.5 is 
generated using P^ and will be called s-pruning front, and the right is obtained using Pj, and 
will be called /-pruning front. In the division of "P, see fig. 8, the non-existing cells below F 
correspond to squares below the /-pruning front, and forbidden cells above F (which occur first 
for a division of V into 2^^ cells) correspond to squares below the s-pruning front. 

As a result of the above construction the two pruning fronts are related. Denote by 8^ and 8f 
the respective coordinates of the two pruning fronts for the same 7. We then find 

^, = i+2(i-^,). (46) 

Moreover the /-front (and the pruned region it encloses) is invariant under Xi, i.e. invariant 
under the map ((5, 7) 1— > (7/2 + 1/2,2(5 — 1). These symmetries in the front induce relations 
between the pruned words listed in table 1. 
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Figure 10: Magnification of the symbol plane around the primary pruning region. The 
s-pruning front left of 6 = 0.5 and the /-pruning front to the right separate allowed 
and forbidden code sequences. The squares of length 2~^, and 2"'' correspond to 
forbidden words of length 10, 12 and 14, respectively. From this one can construct the 
words in table 1. The shown part of the symbol plane corresponds to the lower right 
corner of the Poincare section. The primary pruning region extends to (0, 0.375) to the 
left, to (0,0.75) to the right and up to (0.5,0.5) to the top. 



A complete list of the primary forbidden sequences can therefore be read off from fig. 10, where 
also the corresponding squares with length 2~^, 2~^ and 2"'' are shown. The pruned words of 
width 10,. . . ,14 are given in table 1. 

Looking at fig. 9, one observes that the density of points is much higher close to the /-pruning 
front, while it is low close to the s-pruning front (this is similar to the three-disk billiard 
[21]). The reason is that the mapping from the Poincare section to the symbol plane is not area 
preserving, e.g. the small cells just above F shown in fig. 8 and the much larger ones below F are 
represented in the symbol plane by squares of the same size. Considering a generic trajectory 
which fills the Poincare section uniformly there is only a low probability in the symbol plane 
for squares whose corresponding cell in the Poincare section is small. 



4 Periodic Orbits 

Periodic orbits with period n are fixed points of P", i.e. P^^ = ^. Thus the corresponding 
symbol sequence is also periodic and will be denoted by U, where is a word consisting of n 
letters coi. The same periodic orbit can be denoted by n different words all obtained from one 
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14 


A^BA^B.A^ 









Table 1: Pruned words in the cardioid billiard up to length 14. The shown words 

a.ijj and additionally cu.a, a.Co and a.u are forbidden. Code words which are already 
pruned by shorter words are omitted. Notice, that for example the squares B^A.A^ 
and B^A.A^B of length 12 shown in fig. 10 rule out all words of the form B^A.A^ of 
length 11. 



word by a cyclic permutation of the letters, e.g. ABB = BBA = BAB. Therefore with every 
periodic orbit a class of code words is associated. Two words ZJ and ZJ' belong to the same 
cyclic class if they can be obtained from each other by a cyclic permutation. Periodic orbits, 
which are not multiple traversals of an underlying shorter periodic orbit are called primitive 
periodic orbits. 

For a periodic symbol sequence ZJ = coi . . .cOm the corresponding point in the symbol plane 
((5, 7) is given by 
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where N[ijj) = XlILi ^«'2" ' is the binary interpretation of io. We choose the symbol sequence 
with the smallest value 7 as a representative for that orbit. 

The classical properties length, stability and Maslov index of periodic orbits belonging to dif- 
ferent cyclic classes, which are related by X, T or XT, are identical. Therefore different cyclic 
classes will be combined into symmetry classes if their code words can be transformed into each 
other by refiection and/or time-reversal. Thus it is sufficient to calculate only one periodic 
orbit for every symmetry class. As symmetry classes we define 
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(49) 



using the operations defined in eq. (45) (see [14, 23] for a similar classification in the case of the 
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hyperbola billiard). Within every symmetry class we choose the sequences with the smallest 
binary interpretation N[ijj) as a representative for that class. 

To every symmetry class a multiplicity is associated, which counts the corresponding number 
of different cyclic classes. For Cno the multiplicity is 4, Ct, Cx and Cjx have multiplicity 2 
and Cfuu has multiplicity 1. For each symmetry class at least one example of a corresponding 
periodic orbits is shown in fig. 11. 

In the case of no symmetry ^, ^ and "give the symmetric partners, so that each of the symmetry 
classes is closed under these operations. If U G Cj or U G Cx the corresponding symmetric 
partners are generated by ^ or ^ respectively. For ZJ G Cjx one obtains the symmetric partner 

by ^ or ^, since co =u. Thus despite the fact that the orbit shown in fig. 11c) looks symmetric. 




Figure 11: Fxamples of periodic orbits: a) to = AABABBAB G Cj^y, h) to = 
AAABBB G Cj^u c) = AAABAABBBABB G Ctx, d) to = AABABB G Ct, 
e) = AABAB G Cjf , i) to = AAAABAAB G C|>, g) to = AAAAABAB G Cj| and 
h) cJ = AAABABB G Cno- 
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application ol X yields the corresponding time-reversed orbit. 

A finer classification of orbits with X and full symmetry is obtained by using the symmetry 
lines (see [20] for a similar discussion): 

• C Cx contains periodic orbits of period 2n, which are intersections of the type 
A!2n n . The orbit has two different points with ^ = 0. 

• C Cx contains periodic orbits of period 2n — 1, which are intersections of the type 
n The orbit has one perpendicular intersection with the x-axis and one point 

with ^ = 0. 

• C]| C Cx contains periodic orbits of period 2n, which are intersections of the type X\,^j^^ fl 
x\. The orbit has two perpendicular intersections with the x-axis. 

• ^fuii ^fuU contains periodic orbits of period 4n corresponding to intersections of the 
type n Tq. The orbit has two perpendicular reflections with 90 and one reflection at 
^ = which is traversed in both directions. 

• C\^ C Cfuu contains periodic orbits of period 4n + 2 corresponding to intersections of the 
type X2n^i nTo- The orbit has two perpendicular reflections with dil and one perpendic- 
ular intersection with the which is traversed in both directions. 

Orbits with T or TX or full symmetry are only possible for even n, since " exchanges the 
numbers of A and B. If n is odd, only orbits from C^ and C^o exist. Elements from C^o only 
exist for n > 7. 

In table 2 the number of primitive cyclic classes, primitive symmetry classes and the corre- 
sponding number of elements in the different symmetry classes are shown for n = 2, . . . , 15. 
For period one there are two periodic sequences A and i?, which correspond to the lines J-'± in 
our system. The shortest hyperbolic periodic orbit has period 2 and the sequence AB. 

Due to the pruning not every orbit is allowed. In table 2 the number of pruned symmetry 
classes are indicated by —m for the corresponding symmetry classes. Starting at period 5 the 
orbits A^ B are pruned for j > 4. All other periodic code words up to period 9 correspond to 




Figure 12: Examples of pruned periodic orbits: a) /-pruning: A^B^ b) s-pruning: 
A^B^ and c) combined pruning: A'^BA'^B'^. Note that in c) the reflection in the 
singularity takes place only from above. 
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Table 2: Number of primitive cyclic classes, primitive symmetry classes and the number 
of elements in the corresponding symmetry classes for n = 2, . . . , 15 in the case of a 
complete binary symbolic dynamics. The number m of pruned orbits in the cardioid 
billiard is indicated by — m in the corresponding classes. 



physical orbits. From period 9 the orbits A^BA^B^ i ^ 4 are missing, and for period 10 in 
addition A^BAAB. Three examples of forbidden orbits are shown in fig. 12. 

Some of the forbidden orbits are obviously ruled out by the pruned symbols given in table 1. 
For example the pruned word A^B.A^ rules out the existence of all the periodic orbits of the 
form A^ B, for j > 5 but not for j = 4. The existence of orbits can not be deduced from the 
pruned words (of finite length), but only the nonexistence of periodic orbits. The reason is 
that a periodic orbit is an infinite sequence, which might be pruned by an extremely long word 
which is much longer than the period. This is nicely illustrated in the case of the periodic orbit 
which is forbidden by the word ABA^B.A^BA or BA^B.A^BA^ of length 13, see fig. 10 
and table 1. 

The orbit A^B shown in fig. 12a) is an example of /-pruning. In fact, the corresponding point 
in the symbol plane is inside the pruned square A'^B.A^, see fig. 10. As an example for s- 
pruning in fig. 12b) the forbidden orbit A^B^ is shown, which is pruned by the word B^A.A^^ 
whose square is below the s-pruning front. In addition one also has orbits where combinations 
of both pruning mechanisms occur, see fig. 12c). 

For long orbits especially the s-pruning is the reason for some numerical difficulties in deciding 
whether an orbit has a refiection next to the cusp or whether it is pruned. Some periodic orbits 
close to the cusp have quite large eigenvalues, e.g. for period 11 the most unstable orbit has 
eigenvalue 663393 and the eigenvalue for the most unstable orbit of period 15 is even larger 
than one million. Note that such large instabilities cannot be calculated accurately using the 
monodromy matrix because of round-off errors; therefore we use the method described in [17]. 



20 



4.1 Desymmetrized cardioid billiard 



The code word in the desymmetrized system can be easily read off Irom the code word in the 
lull system in the following way: if two adjacent symbols are identical, assign 0, otherwise 1. 
For example the code for AAAB in the desymmetrized systems is 0011. Periodic orbits of 
the full system with full or TX symmetry are reduced in length by a factor two. For example 
AAABBB G Cmi maps to 001001 = OOT, or AAABAABBBABB G Ctx in the desymmetrized 
system is given by 001101001101 = 001101. 




S 



Figure 13: The rule for the symbols and 1 in the desymmetrized billiard in configu- 
ration space: Connect the cusp with the current point on the boundary. Then refiect 
this line at the normal vector in the current point. If the actual velocity vector is in 
the sector formed by these two lines assign the symbol 1, otherwise 0. 

Likewise one can construct the code word in the full system from the one in the desymmetrized 
system. Only if the code word in the desymmetrized system has an odd number of I's, one has 
to double the word, because otherwise one would not obtain a periodic word. These orbits are 
just the ones with Cfuu or Ctx symmetry in the full system. 

The multiplicities of orbits in the desymmetrized system can be determined from the symmetry 
of the corresponding orbit in the full system: For an orbit from CfuU the multiplicity is still 
1; in the case of Cx and Ct the multiplicity reduces from 2 to 1. Orbits from Ctx and Cno 
are almost indistiguishable in the desymmetrized system, they both have multiplicity two. For 
the relation between the code length n in the full and n' in the desymmetrized system and 
the corresponding multiplicities m and m' the rule ^, ™, = 2 holds. If n/n' = 2 then also 
the geometric length and the stability exponent (which is the natural logarithm of the larger 
eigenvalue of the monodromy matrix M) are divided by two. 

The sign of the trace of M changes with every refiection at the x-axis. The total number of 
refiections at the x-axis is equal to the number of I's in the code word of the desymmetrized 
system. To understand this fact it is helpful to interpret the symbols and 1 in the reduced 
system. 

In the reduced Poincare section the symbol corresponds to all points with \p\ > s/4, and 1 
to all points with \p\ < s/4 (again excluding S and JF). Notice that the region corresponding 
to is not simply connected. In configuration space this translates to the rule illustraed in 
fig. 13). For s < ^16/3 all orbits from the sector 1 immediately refiect at the x-axis. For 

s > ^16/3 only part of these orbits hits the x-axis. The other part with negative p encounters 
a refiection with the boundary dft which gets it closer to s = 0. All the remaining refiections of 
this part get the symbol 0, i.e. they must eventually approach s = and hit the x-axis before 
the next symbol 1 can occur. 
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Table 3: Form of the code words for periodic orbits in the desymmetrized system. 
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Table 4: Form of the code words for symmetric periodic orbits of the full cardioid 
billiard. Y and Z stand for an arbitrary single letter, a for an arbitrary word. 



The fact that the number 7^1 (cu) of I's in the code word gives the number of reflections with 
the X-axis is important for calculating the Maslov index which enters into the Gutzwiller trace 
formula. For a billiard system the Maslov index p is twice the number of reflections at the 
boundary with Dirichlet boundary condition plus the maximal number fi of conjugate points 
[24]. For the desymmetrized cardioid billiard we have in the case of Dirichlet boundary con- 
ditions p = 2n -\- n -\- and in the case of Neumann boundary conditions on the x-axis 
p = 2n -\- n. 

The relation between the code for the desymmetrized billiard and the symmetry classes is easily 
obtained using the symmetry lines. Therefore we first consider the relation between a symbol 
string and the symmetry lines , and Tq in the desymmetrized system [20] . The line 
corresponds to all points in "P, whose symbol sequence is symmetric with respect to the dot 

4 - {^-h ■ (51) 

The line corresponds to all points in "P, whose symbol sequence is symmetric with respect 
to the letter 

A'> ~ {/3.0/3} . (52) 

The line Tq corresponds to all points in "P, whose symbol sequence is symmetric with respect 
to the letter 1 ^ 

To~{/3.1/3} . (53) 

Iterating the basic symmetry lines , x\ and Tq shifts the point of symmetry to the left or 
right. Thus we obtain for the structure of the code words of periodic orbits with a certain 
symmetry the list shown in table 3. 

Now the structure of the code words for a given symmetry in the full system can be determined 
from the structure of the code words for the desymmetrized system. The result is summarized 
in table 4. Notice that if one constructs code words with arbitrary a of the form given for 
Cx, Ctx or Ct, it is possible to obtain a code word with higher symmetry. Furthermore it is 
possible to obtain a code word which is not primitive. The only symmetry class which is not 
incorporated into the above scheme using the symmetry lines are orbits from Ctx- The reason 
is that their code in the desymmetrized system does not have any of the symmetries, like orbits 
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Table 5: Code words of periodic orbits and their length, stability exponent and sym- 
metry class in the full system up to code length 8. 

from Cno- The only difference between them is that orbits from Cno have an even number of 
1, whereas orbits from Ctx have an odd number of I's in the desymmetrized code. Thus for a 
given code length orbits from Ctx in the desymmetrized system occur much more often than 
their number in table 2 for the full billiard might suggest. Actually, in the reduced system 
orbits from Ctx are even more frequent than those from Cno- 

In table 5 we list all periodic orbits up to code length 8 together with their length, stability ex- 
ponent and symmetry class of the full system and the corresponding code in the desymmetrized 
system. 
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4.2 Families of periodic orbits 



Of special interest for the application of Gutzwiller's periodic orbit theory [1] are families of 
periodic orbits which accumulate in length. The code words of the simplest type of families 
have the form A^cu, with u starting and ending in B. One family corresponds to fixed u and 
varying n. Some examples are shown in fig. 14. The family of orbits A^BB yields the shortest 
accumulation length / = 12. Moreover this is the only family of the form A^u with finite cu, 
which exists for arbitrary n. All the other families of this kind are pruned for a finite n > nmax- 

This is nicely seen using the representation of the families in the symbol plane. One has 

2"+™ — 1 2"+™ — 1 
2"+™ — 1 2"+™ — 1 

and thus lim„^oo('^ra, 7™) = ('^^O). with 6 = 2~'^N{u). Furthermore, it turns out, that all the 
points ((5„,7„) lie on straight lines 



N{io) 
In = — 7^ 



{2"^S^-N{u)) . (56) 




Figure 14: Fxamples of members of some families of periodic orbits which accu- 
mulate in length a) A'^BB for n = 11, b) A^BBB for n = 11, c) A^BABB for 
n = 11, d) A^BA B^AB for n = m = 11, e) A^BABAB^AB for n = 7, m = 6 f ) 
A^B^A'^BAB^AB. 
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Figure 15: Sequence of points [6^ 7^) in the symbol plane of some pruned accumulating 
families A'^co together with the pruning front. The labels denote the corresponding 
shifted code word of the first plotted member. 



By shifting the code word or by considering a symmetric partner it is possible to obtain 6 G 
]0.375, 0.75]. If different words fulfill this criterion the one with the smallest distance to 6 = 0.5 
is chosen. 

In fig. 15 some examples are shown together with the pruning front. From the figure one can 
read off the last allowed orbit as well as the pruning mechanism. If 6 g]0.375,0.5[ the family 
is s-pruned, whereas for 6 g]0.5,0.75[ the family is f-pruned. As explained above it is also 
possible that a family is pruned by both mechanisms. 

We observe that the pruning front is a monotonously increasing curve in the interval [0.375, 0.5], 
and monotonously decreasing for [0.5,0.75]. Moreover it only reaches at 0.375 and 0.75. 
Therefore families of any kind can only exist if at least one of their limit points in the symbol 
plane is (0.75,0). In particular we think that the family A^BB exists for any n. 

One can construct a variety of more complicated families, e.g. with the structure or 
see fig. 14. The special case A'^BAB'^AB^ see fig. 14e), has the accumulating length 
/ = 24 and probably exists for any n, because the corresponding limit point in the symbol 
plane is (0.75,0). Actually we think that there is an infinite number of accumulating families. 
They can be constructed by introducing two combined symbols Fn = AB^A and Gm = BA^B. 
Now any sequence of F^s and Gm's gives an accumulating family, as long as all the indices 
collectively go to infinity. One example with different rii is shown in fig. 14f). The length of the 
limit orbit of such a family is given by the total number of the symbols Fn and Gm multiplied 
by 12. 

The stability of periodic orbits from accumulating families grows with increasing period, because 
they must come arbitrary close to the cusp. However, their stability grows much slower than 
for generic orbits because they also get arbitrary close to the parabolic line JF. 



25 



The pruning of the "almost accumulating" families is related to the limiting orbit which almost 
hits the cusp. It consists of a sliding motion and of a finite orbit hitting the cusp as illustrated 
in fig. 14. We now turn to the study of these finite orbits. 

4.3 Finite orbits 

In the quantum mechanical billiard problem discontinuities of derivatives of the boundary curve 
can play an important role (see, e.g. [25, 26, 27] and [28, 29, 30, 12, 31] and references therein). 
The cusp orbits start at an arbitrary angle in the cusp and eventually return to the cusp with 
an arbitrary angle without the need to fulfill any refiection condition in the cusp. We think 
that the investigation of orbits having a refiection in the cusp is already important from the 
classical point of view: on the one hand due to the connection to families of periodic orbits 
and on the other hand because cusp orbits are also extrema of the Lagrangian, as it is the 
case for periodic orbits. This can be seen in the following way: the variation of the Lagrangian 
C = C{cf>\, . . . , cf>n) yields a system of n coupled nonlinear equations V£ = 0. If cf>i = zLtt the 
corresponding equation is automatically fulfilled. 



regardless of the values of cf>i-\ and ^i+i. This means that there is no refiection condition to 
fulfill in the point cf)i = zLtt. Thus in some sense the singularity looks like an infinitely small 
circle, i.e. it can refiect into any direction. Notice that the above argument only holds if we 
look at £ as a function of cf)i. This dependency of the Lagrangian on the choice of coordinates 
is due to the singularity. We will not consider cusp orbits as periodic orbits, but instead as 
finite orbits with the same initial and final point in configuration space. 

In order to classify the cusp orbits we extend the symbolic dynamics with a third letter C, 
which corresponds to the lines S and V . For a cusp orbit starting in S the next point lies on 
F~^, from which one iterates until the first point lies on F. Since the momentum is not defined 
in the next step we consider cusp orbits as finite orbits whose code word is CaC, where a is a 
code word consisting of A and B. One could omit the C's, but we find it more convenient to 
have the number of letters in a code word correspond to the number of refiections. Thus the 
shortest cusp orbit, which runs along the symmetry line has the code CC, because it starts in 
S and has its last refiection on F, see fig. 16a). 

Since the partition of the Poincare section as shown in fig. 6 and fig. 7 is given by the iterates 
of F, cusp orbits correspond to the intersections of these lines. Actually it is simplest just to 
iterate F~^ and to look of intersections with F. After n iterations we have 2""'"^ — 1 intersections 
(ignoring pruning) i.e. at every iteration there are 2" new intersections, which is exactly the 
number of cusp orbits one expects without pruning. Defining F~"~^ = p^P"^ the number of 
refiections of a cusp orbit which is given by the intersection of F~"~^ and F, is n + 2. 

Similar to the case of periodic orbits one can classify finite orbits hitting the cusp according 
to their symmetry by using the symmetry lines. Intersections of Tq and F~" correspond to 
T-symmetric cusp orbits which have period 2n because the distance to F and F~^ is n — 1 
forward respectively backward iterations. Intersections of and F~" yield X-symmetric cusp 
orbits. In the same way cusp orbits with Xi symmetry are found. If one starts on a symmetry 
line past and future are symmetric and thus the symmetry is reproduced in the cusp, i.e. a 



dC_ 








(57) 
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Figure 16: Examples ol finite orbits a) CC, b) C AC , c) C AAC , d) C ABC , e) C AAC , 
1) CAABC, g) CAB AC, h) CAAAAC and i) 



T-symmetric orbit retraces itsell also in the cusp and a X-symmetric orbit has a X-symmetric 
refiection in the cusp. Note that the latter use the vertical tangent in iS, while the others use 
any tangent vector. The exceptional period two cusp orbit CC has full symmetry. This is not 
the case for any other cusp orbit. 

By counting the number of intersections of the symmetry lines and iterates of V (ignoring 
pruning) we find that there are 2("~^)/^ symmetric cusp orbits of even period n and 2("~"^)/^ for 
odd n. Like in the case of periodic orbits the non-symmetric cusp orbits are the overwhelming 
majority for large n. For a complete symbolic dynamics the number of cusp orbits is exactly 
2"~^ . where n is the total length of the word C aC . So already this number grows faster then the 
number of periodic orbits because we do not have to consider cyclic classes here. Furthermore 
cusp orbits can be combined to give multiple cusp orbits. Their number can be estimated 
using the topological polynomial [32] for the symbolic dynamics with letters {A, 5,(7 = CC}, 
which is \ — — , such that we obtain /itop = ln(l + a/2) for the topological entropy. Thus 
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they occur much more often than the periodic orbits for which the htop = In 2 in the case 
of a complete symbohc dynamics. Actually we do not gain any information by this because 
every multiple cusp orbit can be uniquely decomposed into single cusp orbit, but classically 
they become interesting if we smoothen the singularity and quantum mechanically they are 
interesting as multiple diffractive orbits. Finally note that in considering multiple cusp orbits 
one can produce full symmetry and TX-symmetric cusp orbit by superposition of single cusp 
orbits with their appropriate symmetric partners. In this way we can even build multiple cusp 
orbits with any symmetry starting with a cusp orbit with no symmetry. 

We will now briefly explain the relation between cusp orbits and the families of accumulating 
orbits. The families A'^co perform a large number of successive reflections along the boundary, 
which are followed by some other reflections after which the orbit returns to the approximate 
sliding motion. For large n this family tends to the sliding motion plus a limiting cusp orbit. 
If A^BaB is the code word for the family, the corresponding cusp orbit is given by CaC . The 
largest angle x between the x-axis and the limiting cusp orbit determines the maximum number 
of consecutive reflections having the same letter. Only for x = a family can exist for any n. 
Otherwise we have a finite family of short orbits which do not accumulate. 



5 Global properties of the billiard 

In this section some global characteristics of the cardioid billiard are calculated. In the first part 
we obtain an analytical result for the average distance between refiections / and an analytical 
estimate for the Kolmogorov- Sinai entropy /^ks- The same quantities are then calculated using 
averages over periodic orbits, where each orbit is weighted by its stability. 

We first calculate dynamical averages by iterating the billiard map. Averaging over several 
generic trajectories for 5 X 10^ iterations we obtain for the average length between refiections / 
and the Lyapunov exponent, which in our case approximates the KS entropy /^ks 

/ ^ 1. 851 (58) 
/iKS ~ 0.653. (59) 

/iKS is calculated according to [33], however, we use the linearized map (16). One can also 
calculate these quantities analytically: 

For any billiard in the domain the average length is determined by the area |0| and circum- 
ference |90| by transforming the average over the length function /(s,p) on V into the volume 
of the energy shell [33] 

/ = J^l{s,p)d^ = ^^^^ J^l{s,p) ds dp (60) 
= , , / dx dy dS = | ^1 , , i.e. for the cardioid (61) 

2|ao| ifix5i ^ ^ \dny ^ ' 

= -^TT^ = 1.85055... (62) 
16 

where dfi = ^j^dsdp is the normalized Liouville measure on V. 

An analytical expression for a lower bound on h^s was given by Wojtkowski in [3]: 

hKS > / log d/i , (63) 
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where q estimates the contribution to the expansion from below. Since q contains the length 
l{s,p), which is given by the root of a cubic equation we again pass to the generating function 
of the map /(^i, ^2)- The variables are transformed by 

5((^i,(^2) = 4sin((^i/2) (64) 

(A A\ I (iv(^l,^2),r(^l)) ir... 

P{<pi,<p2) = {v,Ti) = . (65) 

Since s is independent of (f>2 we obtain for the Jacobian of the transformation 

d{s,p) ( {L,T){L,L') \ 

where a prime denotes the derivative with respect to ^2- Since not all combinations (^1,^2) 
correspond to allowed trajectories, we have to restrict the integration ranges. Moreover we use 
the symmetry of the system to restrict the integration to positive ^1 and obtain 

hKS>T^/[ J gi^r,h)det dM^,. (67) 

\aii\J J 0(^1,^2) 

U — 7r+0i 

Let the linearized Poincare map be given by 

^^=f! (68) 



with det DP = 1. The arguments of Wojtkowski [3] for matrices with only positive (or negative) 
entries work equally well for checker board matrices (18). Then g = \fad \ \fhc. Using (16) we 
obtain = {k\ — 1)(A;2 — 1) where hi = tui/rii. Therefore we have to evaluate the following 
integral 

hKS> J J log (^^/{h - l){k2 - 1) + \/hk2 -h- k2^ det ^^^^d(^2d(^i . (69) 

— 7r+(/)i 

The numerical integration gives h^s ^ 0.633. Even though Wojtkowski's theorem only gives a 
lower bound on h^s this value is surprisingly close to the numerically measured value. 

Uks can also be approximated by calculating the entropy of the partition of the Poincare section 
V into cells labeled by words of length n 

hKs{n) = --J^P'^^^P'^^ i^^) 

where = /u(cu) is the size of the cell with label u [34]. Since the cell sizes in V are quite hard 
to obtain we determine them numerically in the symbol plane. This is done by calculating the 
probability with which each cell is visited by performing many iterations of one initial condition. 
Thus we numerically approximate the invariant measure in the symbol plane by assigning a 
probability to every cell labeled by a word u of length n. Thus the sum (70) extends over 
all allowed cells in the symbol plane. If we denote the total number of allowed cells labeled 
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by words of length n by N[n) = Y,uj 1? we obtain an estimate for the topological entropy with 
respect to code length with p = 1/N[n) 

/itop(f^) = --V'plnp = -lniV(n) , (71) 

For n = 22 we obtain h^s ~ 0.66 and htop ~ ln(1.98) ~ 0.68. h^s is smaller than /itop! 
the difference tells us how much the invariant measure on the symbol plane deviates from 
equipartition. The three values for h^s are 0.63 < 0.65 < 0.66. Their ordering is correct 
because the analytical value is a lower bound and the value obtained from the symbol plane 
tends to overestimate due to finite size effects. 



5.1 Global properties from periodic orbits 

Now we repeat the calculation of /^ks, ^top and / using periodic orbits. The main point in 
the calculation of these averages is that the sum in (70) can be transformed into a sum over 
periodic orbits of fixed length. Each cell visited by a periodic orbit 7 is assigned a probability 
given by the inverse of the eigenvalues e""'', from which the invariant measure is approximated 
(see, e.g., [35, 32, 36]). 

This approach can be illustrated by a plot of the distribution of periodic orbits in the Poincare 
section (see fig. 17). We observe that the distribution is not as uniform as one might expect. 
Most obvious we have regions with very few periodic orbits in the neighborhood of JF, because 
these orbits are rather stable. The same structure can be observed in the division of V into 
cells in fig. 7. In a region where the cells are small there is a high density in contrast to regions 
with large cells, where the density is low. As we already discussed the accumulating families are 
not only close to but also close to S and therefore close to F and its images and preimages. 
These families are the only ones that evade the region surrounding F exactly at the points 
where there exists a corresponding cusp orbit. Parts of F with large s have a high density 
of periodic orbits in their neighborhood. They come close to F without having been close to 
JF, and are therefore extremely unstable, which is the reason for their high density. Moreover 
some of the basic symmetry lines and their iterates are visible. The reason for a relatively high 
density of points along the symmetry lines is the fact that they are one-dimensional. Thus the 
small fraction of symmetric orbits has to fill only a "small" set in "P, whereas the large number 
of orbits without symmetry has to fill large areas in V. Plotting only orbits of higher period 
the symmetry lines are less visible. If one plots the periodic orbits in the symbol plane, they 
are much more uniformly distributed, see fig. 18. 

The weighted average over periodic orbits of length n of a quantity quantity /(u, /) is now given 
by 

if): = ^^{1""^?""^ , (72) 

where the sum runs over all periodic orbits 7 of code length n, including multiple traversals. 
Although for large n we do have E7 e""'' ~ 1 [37] we include this term for normalization. 

Using the periodic orbits up to code length 20 we we calculated the mean stability exponent 
per refiection {u)'^/n. The values approach a constant for increasing n with (m)^q/20 = 0.657, 
which agrees quite well with the the value of h^s in (59). 
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Figure 17: Plot of the points of all periodic orbits with code length up to n = 15 in 
the Poincare section V . The symmetry lines and the are visible, comp. to fig. 4. 

For the average distance between consecutive refiections we obtain (/)^q/20 = 1.864, where 
again the values are slowly decreasing for increasing n, such that we have good agreement with 
the analytical result (62). In fact the conversion of values measured with respect to the discrete 
time (Poincare section) and the one measured in phase space is just given by the factor (62). 
Thus the corresponding quantities for the continuous fiow are given by 

/iKs(20)// = 0.36 (73) 

^{urjl = 0.36. (74) 

Of course h^s/l can also be approximated by averaging u/l which gives (m//)^q = 0.363. 

For the family of billiards (1) Robnik calculated the KS-entropy in [7]. He gave values for e 
close to 1, from which we extrapolate h^s = 0.34 for the cardioid, which is consistent with our 
values. 

Finally we calculate the topological entropy htop from the growth of the number of periodic 
orbits by considering htopin) = MnA^(n), where N{n) is the total number of points belonging 
to all the periodic orbits of period n, including multiple traversals. We obtain a plateau for 
n > 10 with htop ~ 0.683 ~ In 1.98, which is in good agreement with the value obtained in the 
symbol plane. 

5.2 Statistical properties of periodic orbits 

The following investigations are based on periodic orbits up to code length 20 (See, e.g. [14, 
23, 38, 39, 40] for similar studies for other systems). The principle difference of the averages 



31 



0° taken in this section to the averages O"" calculated in the previous section is that now we 
take all periodic orbits as equally weighted, i.e. we have a uniform weight in the symbol plane. 
Thus we transfer the picture of periodic orbits in the Poincare section (fig. 17) into the symbol 
plane (fig. 18) and average /(u, /) over periodic orbits 7 of code length n by 

(/): = ]^E/K,M- (75) 

Looking at fig. 18 we first observe that in fact the orbits are distributed rather uniformly in the 
non pruned region. The symmetry lines are almost invisible. Only the accumulating families 
show up as holes that are not related to the pruning front. 




Figure 18: Plot of the points of all periodic orbits up to code length n = 15 in the 
symbol plane, compare to fig. 17 where the corresponding points are presented in the 
Poincare section. 

In the last section we calculated the growth behaviour in terms of the code length, but also 
the growth behaviour of the number Af{l) of periodic orbits 7 with geometric length below 
a given length / is of interest. Because the cardioid billiard is strongly chaotic one expects the 
typical exponential proliferation 

A/'(/) ~ Ei (r/) , (76) 

Tl 

where r is the topological entropy with respect to the geometric length of orbits. However, this 
spectral staircase is not well defined in our case because of the families accumulating in length 
(the same happens, e.g. in the case of the stadium billiard or the wedge billiard [39]). 

Therefore we define Sacc as the set of all periodic orbits 7 whose corresponding code has more 
than five consecutive letters A or i?, and Sj-gg as the set of periodic orbits with up to five 
consecutive A oi B. This choice is somewhat arbitrary, but motivated by the observation that 



32 



orbits from Sacc already have the geometrical structure in configuration space of the possible 
limit orbit of the considered family. We define A/'reg(0 as the counting function of the number 
of orbits from Sj-gg with geometric length less than /. In fig. 19 A/'(/) and Afj-eg{l) are shown 
in logarithmic representation using the periodic orbits up to code length 20 together with a 
fit with the asymptotic behaviour (76) r 0.345 for A/'reg(0- One clearly observes for Af{l) 
the strong increase next to length / = 12, which is caused by the family A'^BB. This step 
would even be more pronounced, if orbits of higher code length had been used. In contrast the 
logarithmic plot of A/'reg(0 is a much "smoother" curve, with most of these "steps" removed. 
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Figure 19: In a) the logarithmic representation of the number of periodic orbits with 
geometric length less than / is shown using all orbits up to code length n = 20. In 
b) A/'reg(0 is shown in logarithmic representation where all orbits with more than five 

tI 

consecutive A or B are excluded; also shown is ^ for r = 0.345. 



That the above splitting allows one to select orbits belonging to accumulating families is also 
visible if one plots the distribution Pn{l) of lengths for a given code length. In fig. 20a), where 
all orbits with code length 14 were used, long tails down to / = 10 are visible. In contrast in 
fig. 20b) only orbits from Sj-gg were used, and the outliers are removed. 

We calculated the average length of periodic orbits and observed the expected increase of the 
mean length (/)° = nl with the code length n, with / = 1.99 using all orbits and / = 2.054 
using orbits from Sj-gg- The variance of the distribution increases with n. For orbits from ^reg 
one finds that cr^ a^n for n > 15, with 0.47. Furthermore we observed, as in the 

case of the hyperbola billiard [14], that Pn{l) are approximately Gaussian distributed. This is 
demonstrated in fig. 20d) where each Pn{l) is shifted to the origin and the variance is normalized. 

In [14] a relation between the number of primitive cyclic classes of the code words and the 
classical staircase was derived under the assumption that Pn{l) has a Gaussian distribution and 
that the mean length (/)° and the variance cr^ depend linear on n. The result adopted to our 
case is 



a2 



h 



top 

T' 



1 + a/1 - 2htop^yp 



(77) 



where htop denotes the topological entropy with respect to code words. In the case of orbits 
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Figure 20: Probability distribution Pn{l) for n = 14: a) all orbits are included, b) only 
orbits with not more than five consecutive A oi B are used, and in c) the difference 
between both is shown. In addition in a) and b) a Gaussian distribution is shown 
as dashed line. In d) the shifted and normalized distributions Pn{l) are shown for 
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,20. 



from ^reg we obtain htop = 0.672 and therefore 

T ^ 0.34 , (78) 

which is in good agreement with the value r 0.35 obtained from the fit of Afj-eg{l) with the 
asymptotic behaviour. The consistency of these results can be seen as justification of the above 
splitting. 

We also considered the distribution of the stability exponents for fixed code length n and 
found that is centered around a mean = nu with u = 0.70 using either all orbits or 
only orbits from Sj-gg- The distributions are stronger peaked than Gaussian distributions. 

A further common statistics is the spacing between neighbouring lengths. In order to obtain 
a mean spacing of one, the asymptotic behaviour is used to unfold the length spectrum. The 
result using orbits from Sj-gg with length / < 30 is shown in fig. 21 is in agreement with the 
Poisson distribution P{s) = e~*. 

It is also interesting to look at the dependence of the stability exponents vs. the geometric 
length l^. The result is shown in fig. 22. Some of the families are clearly visible as series of 
points with accumulating and thereby increasing u^. 
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Figure 21: Length spacing distribution P{s) (full line) together with the Poissonian 
distribution e~* (dashed line). 
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Figure 22: Plot of the stability exponents vs. the geometric length using all orbits 
up to code length 20. 

6 Summary 

In this paper we have developed a symbolic dynamics for the cardioid billiard by constructing 
a partition of the Poincare section V using two symbols A and B. After completion of the 
part on the symbolic dynamics a preprint appeared [12], where the same coding was found 
independently and used to study diffraction effects of the quantum mechanical system. 

For the symbolic dynamics it turned out that not every sequence of symbols is allowed. We 
obtained the pruning front, which comes in two pieces related to the two pruning mechanisms 
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in the system. Using the symbohc dynamics periodic orbits can be labeled in a unique way. As- 
sisted by the knowledge that all orbits correspond to maxima of the Lagrangian we calculated a 
large number of periodic orbits up to code length 20. Complete sets of higher periodic orbits are 
hard to obtain because some of them are extremely unstable due to the unbounded curvature. 
Using the symmetry lines of the billiard map, we presented a classification of periodic orbits 
with respect to their symmetry properties. This was first obtained in the symmetry reduced 
symbolic dynamics and then translated back to the full system. Combining Wojtkowski's result 
about convex scattering billiards with a geometric argument in the symmetry reduced system 
we were able to determine the Maslov indices from the code words. 

Studying families of periodic orbits with short geometric length provides a nice application 
of the pruning front, because it allowed us to determine whether a family exists for arbitrary 
code length, or whether it is eventually pruned. In the latter case a plot of the points of the 
periodic orbits in the symbol plane allows for the determination of the last allowed member of 
the family. The converse argument enables us to write down an infinite number of families that 
accumulate in length. Furthermore it turned out that finite respectively cusp orbits appear 
rather natural even from a merely classical point of view as parts of the possible limit orbits of 
the short families. 

In the last section we calculated an estimate for the Kolmogorov- Sinai entropy and found good 
agreement with numerically calculated values. Averaging the periodic orbits we find consistent 
values for the KS entropy and the average length between refiections. We illustrate the idea of 
the periodic orbit averaging by a plot of periodic orbits in the surface of section. The topological 
entropy is quite close to In 2 because the pruning sets in rather late. 

To obtain well defined statistics despite the presence of accumulating families we suggested 
a method to subtract the accumulating families. This procedure might also be helpful when 
using Gutzwiller's periodic orbit theory [1], where families of accumulating periodic orbits have 
to be treated separately. We believe that this method might be useful for other systems with 
accumulating families. 

With the complete set of periodic orbits, knowledge of the Maslov indices and an understanding 
of accumulating families and cusp orbits all the ingredients for the periodic orbit quantization 
of the cardioid billiard are available. Eventually the same program should be carried out for 
Robnik's family of billiards, at least close to the cardioid, where it is conjectured to be ergodic. 
On the one hand, the system then lacks the singularity which makes it simpler on first sight. On 
the other hand we think that it becomes more difficult, because there will be inverse hyperbolic 
orbits, such that it is not sufficient to look for maxima of the Lagrangian in order to find 
periodic orbits and moreover a binary symbolic dynamics will not suffice. 
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